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Three static models with two interacting phantom and ghost scalar fields were considered: a 
model of a traversable wormhole, a brane-like model and a spherically symmetric problem. It was 
shown numerically that regular solutions exist for all three cases. 
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^ ' I. INTRODUCTION 

(N 

\ Scalar fields play a fundamental role in the modern cosmology and astrophysics. Having come in them from the 
• theory of elementary particles, the scalar fields are widely used both at creation of models of compact objects and 
P^h ' at research of the evolution of the Universe as a whole. In particular, the models of quasi-star objects from scalar 
fields, the so-called boson stars, are well-known [H, [2>[3[. In turn in cosmology scalar fields are the basis for creation 
of models of the early inflationary Universe 

Another field of application of scalar fields is connected with one of the most excitingevent in astrophysics and 
cosmology of the last decades - discovery of the acceleration of the present Universe [H, Q . This discovery has 
stimulated an appearance of a large amount of various models trying to explain this phenomenon. A basis of all the 
models consists in violation of different energy conditions. For cosmology, most essential are two energy conditions 
(see, e.g., Ref. [(|): the first one, the so-called strong energy condition (SEC), states that p + 3p > 0, where p and p 
are effective energy density and pressure of matter which determine the evolution of the Universe. In hydrodynamical 
language this means that the parameter of the equation of state w > —1/3. This case corresponds to the decelerated 
expansion of the Universe (the Friedmann models). Violation of the SEC leads to the accelerated expansion of the 
Universe (the inflationary models). The second condition, the so-called weak energy condition (WEC), states that 
p + p > 0, whence it appears that w > — 1. Failure satisfying the WEC results in exponential or more fast accelerated 
■ expansion of the Universe. The substance providing such acceleration was called dark energy. This term usually means 
that the equation of state lies in the range — 1 < w < — 1/3 (the cosmological constant, quintessence [7|, Chaplygin 
gas Q, the theories of gravitation with high derivatives Q, etc.). All the models describe the present accelerated 
' expansion of the Universe within the accuracy of observations. 

The mentioned above violation of the SEC is normally arisen at some choice of a potential energy of usual scalar 
fields (see examples of such potentials in Refs. 6]). However, as it was shown in Ref. [l0| by a model- independent way 
on the basis of study of data sets containing 172 SNIa, it is possible in the present epoch that w < —1, i.e. the WEC 
can be violated. (For more recent papers with the latest observational data see Refs. [HI].) Such an unusual state 
of matter is also known as phantom dark energy [l2j . There are models of phantom energy within the framework of 
higher-order theories of gravity [l3[ , braneworld cosmology [l4[ , etc. 

Another popular direction consists in consideration of ghost scalar fields with negative sign before kinetic term [l5| . 
In Ref. [l6| we considered the cosmological model of the early Universe with two interacting ghost scalar fields with 
special form of the potential energy: 

X ) = y (0 2 - ml) 2 + y(x 2 - m\f + 0V - V . (1) 

Here cj>, \ are two scalar fields with the masses m\ and m,2, Ai,A2 are the self-coupling constants and Vq - some 
constant which should be chosen on the assumption of a problem statement. The essential feature of this potential is 
presence of two local minima at % = 0, (j> = ± mi- It allows existence of localized solutions with finite energy in 
problems with such potential. As is known [l7|, for a case of one scalar field, localized solutions could exist only for a 
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case with two or more minima of potential energy when solutions start in one minimum and tend to another one. In a 
case of two scalar fields, a situation is possible when a solution starts and finishes in the same minimum. In Ref. fiH ] 
such a case was carried out: the solution started at t — > — oo and returned to the same local minimum at t — > +00. 
Such type of soliton-like solutions, localized on a spacelike hypersurface, known as spacelike branes (S-branes) (l8j . 

The mentioned in Ref. |16j possibility of existence of the localized in time solutions for phantom fields allows us 
to hope for presence of similar solutions and for a static case. In this paper we consider three models created by 
two interacting phantom and ghost scalar fields with the potential (Q]): 1) a traversable wormhole; 2) a brane-like 
solution - analog of a domain wall solution but with asymptotically non-flat spacetime (anti-de Sitter spacetime); 3) 
a spherically symmetric particle-like solution. 

The model of so-called traversable wormhole was suggested in Ref. (for a recent review, see [20]). This wormhole 
is created by some special matter with the violated WEC. As phantom scalar fields also violate the WEC, they could 
be used at modelling of traversable wormholes. The researches of such models were already carried out (see, e.g., 
Refs. [U HH). In these works some effective hydrodynamical energy- momentum tensor with the equation of state 
w < — 1 was chosen as a source of matter. But a distribution of this matter was added by hand and, correspondingly, 
the nonself- consistent models of the traversable wormholes were considered. In this paper we consider a self- consistent 
model of a traversable wormhole created by two interacting phantom and ghost scalar fields with the potential (fl}. 
The specified above features of this potential allow us to find regular solutions with localized energy density. 

Domain walls are topological defects and they arise in different aspects both in particle physics and cosmology (see, 
e -g- ; [4. 123I l24j and references therein). They separate a spacetime into several domains along a single coordinate. In 
a case of scalar fields domain wall solutions exist when the scalar filed potentials have isolated minima. The domain 
walls are surfaces separating minima of the potentials with different vacuum expectation values. The region of fast 
change of the scalar field corresponds to the domain wall. The domain wall refers to as thin wall if the energy density 
of the scalar field is localized at the domain wall surface and could be re plac ed by the delta function. The different 
variants of the domain wall solutions were found in Refs. [IH [26|, [U [H, with asymptotically flat, de Sitter and 
Schwarzschild spacetimes. There are also known the so-called thick domain wall solutions [30, bLU, [32, [33| which could 
exist at late-time phase transitions in the evolution of the Universe. In this paper the consideration of the thick 
domain wall model with the potential ([T]) is presented. 

Search of spherically symmetric solutions with various matter sources always was an important problem in special 
and general relativity. Such solutions are using both at investigation of different particle-like models of elementary 
particles and creation of models of star-like objects and another large-scale configurations. The source of matter is 
fields with various spins both interacting between each other and with gravitational field (see, e.g., Ref. 
are well-known regular solutions for the scalar fields both for noninteracting and self-interacting fields [j], [ 
case of usual (non-phantom) scalar fields, the model of boson star with potential fl]) was considered in Ref. 
was shown there that there are regular solutions in the case under consideration. Recently the spherically symmetric 
model with self-gravitating matter with the hydrodynamic equation of state w < — I was considered in Ref. [35| . 
Below we show that regular solutions for spherically symmetric case exist for phantom and ghost scalar fields also. 

The paper is organized as follows: in section [IT] the general gravitational and field equations for all three above 
mentioned models are presented. In sections IIII1 IIVI and [V] the models of traversable wormhole, brane-like and 
spherically symmetric solution are considered. In section |VT1 we present comments and conclusions. 




II. GENERAL EQUATIONS 



We chose the Lagrangian as follows: 



L = 



R 



16nG 



(2) 



where R is the scalar curvature, G is the Newtonian gravity constant and the constant e = ±1. In the case e = +1 
one has the theory of usual scalar field plus gravitation. The case e = — I corresponds to the theory of ghost scalar 
field. The corresponding energy-momentum tensor will then be: 



T 



and variation of the Lagrangian @ gives the gravitational and field equations in the form: 

G k ; = 8TrGTf, 



(3) 



(4) 
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d(<p,x)' 



(5) 



In our case equations (01)- © are the system of ordinary nonlinear differential equations with the potential energy 
from ([T]) . As it follows from the experience of previous researches of similar systems [16| , finding of solutions of the 
system (HI)-© * s reduced to search of eigenvalues of the parameters mi, m 2 . Procedure of search of solutions and its 
application for investigation of the models mentioned in Introduction will be considered in next three sections. 



III. TRAVERSABLE WORMHOLE 

We will search for static solutions of equations (HI)-© for the following metric [36| : 

ds 2 = B{r)dt 2 - dr 2 - A(r)(d8 2 + sin 2 Odcf 2 ), (6) 

where A(r), B(r) are the even functions depending only on the coordinate r which covers the entire range — oo < r < 
+oo. Using this metric, one can obtain from © and © the following equations (at e = — 1): 



A" 1 /A'\ 2 lA'B' 



A 2 \AJ 2 A B 

at_ i^^_i(^y _l(&\ 3 

~A + 2l"B 2\A) ~2\BJ + ~B 

1 [A'\ 2 1 lA'B' I 



= 2 



<P' 2 + X.'\ (7) 
1 



2 y 2 + x' 2 ) + v 



(8) 



where a prime denotes differentiation with respect to r. Eq. ([7]) was obtained by subtracting (£) component from 
(*) component of the equations ©, and the equations © and © are (g) and (£) components of Eqs. ©. The 
corresponding field equations from © wm be: 

^ + (t + 5¥ ) ^ = v [2x ' + Al( ^ _ ro?)] ' (10) 

^ + ^f)x' = x[2^ + A 2 ( X 2 -m|)] . (11) 



In the equations (|7])-([TPj) the following rescaling are used: r — * yj8irGr, ip — > </?/v87rG, \ ~ * x/v8ttG, toi, 2 — > 
toi^/vottG. 

As it was shown in previous researches of problems with the potential (JTJ) (see, e.g., [3), regular solutions of the 
system of nonlinear differential equations ([7|)-(|lip could exist only for some values of the self-coupling constants Ai, A 2 
and the masses of the scalar fields mi , m 2 , and also depend on boundary conditions which set under the problem 
statement. Particulary, specifying some values of the parameters Ai,A 2 , one has already effect on the shape of the 
potential |T]) that, in turn, determine a possibility of existence of regular solutions of the system (fT])- (fTTj) . The further 
task consists in a search of such parameters mi,m 2 which give regular solutions. In this sense the problem reduces 
to a search of eigenvalues of the parameters mi, m 2 and corresponding eigenf unctions A, B, tp and x f° r the nonlinear 
system of differential equations (|T])- (fTTj) . 

The technique of solution of systems similar to (fT])- (fTTj) is described in Ref. [37j in details. The essence of this 
procedure is the following: on the first step one solve the equation (|10[) with some arbitrary selected function \ looking 
for a regular solution existing only at some value of the parameter mi. On this step the influence of gravitation is 
not taken into account. Then this solution for the function tp insert into the equation (|11|) and one search for a value 
of the parameter m 2 yielding a regular solution. This procedure reiterate several times (three usually enough) for 
obtaining of acceptable convergence of values of the parameters mi,m 2 . The obtained functions tp, \ are inserting 
into the gravitational equations ([T])-©. Eq. ©, which is the constraint equation, is using for specifying of boundary 
conditions (see below). The obtained solutions for the metric functions A, B are inserting into the complete equations 
for the scalar fields (fTU|) - (fTT]) and they are solving again for a search of eigenvalues of the parameters mi,m 2 with 
account of gravitation. This procedure reiterate so many times as it is necessary for obtaining of acceptable convergence 
of values of the parameters mi, m 2 . 

The described procedure of a search of solutions of the system (|T))- (fTTj) . also known as the shooting method, allows 
to find rather fast values of the parameters mi,m 2 at which regular solutions exist. We have checked the obtained 
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FIG. 1: The scalar fields ip, \ in the wormhole model for the 
boundary conditions given in (|12[) . 



FIG. 2: The energy density Tq (r) for the wormhole model. 



solutions using the NDSolve routine from Mathematica substituting the eigenvalues mi, m 2 and solving (|T |) . ([8]) . fllOjl . 
(HU) directly. 

The boundary conditions are choosing with account of 7L 2 symmetry in the following form: 

ip(0) = n/3, 
X (0) = VM, 

1 



A(0) 
B(Q) 



F(0(O),x(O)) : 
1.0, 



<pf(0) = 0, 

x'(o) = o, 

A'(0) = 0, 
B'(0) = 0, 



(121 



where the condition for A(0) is choosing to satisfy the constraint © at r = 0, V(<j>(0),x(0)) is the value of the 
potential at r = and the self-coupling constants Ai =0.1 and A2 = 1. 

Then, using the above procedure for obtaining of solutions of the system ([?])- (1TT1) , we have the results presented 
in Fig. QJ-f!]). These results are obtained for the masses mi « 2.661776085 and m 2 ~ 2.928340304. As one can see 
from Fig. [[Tj). <p — > mi and % — > as r — > ±00. It corresponds to asymptotic transition of the solutions to the local 
minimum of the potential Jl} (see Introduction) . The arbitrary constant Vq was chosen in such a way that the value 
of the potential in the local minimum was equal to zero, viz Vq = (A2/4)m|. Such a choice of Vq ensures zero value 
of the energy density as r — * ±00 (Fig. $Z§). 

Let us estimate an asymptotic behaviour of the solutions. For this purpose we will seek for solutions of the equations 
PU ]) -([TT ]l in the form: 

<p = mx-6tp, x = $X, (13) 
±00. Then the right hand side of Eqs. d)-© goes to zero and their particular solutions 



where Sip, Sx <C 1 as r 
are: 



A 
B 



2 1 2 



B a 



1 



(14) 
(15) 



where ro and B^ are constants. Practically Tq defines a total mass of a wormhole and B^ the run of time at the 
infinity. By corresponding redefinition of the time t, these solutions could be reduced to a flat form in spherical 
coordinates, i.e. one has asymptotically flat Minkowski spacetime (see Figs. ©J®)- Then, taken into account 
(fT4]) -([T5 j) . the corresponding asymptotic equations for the scalar fields (HQj) - (jTTJ) will be rewritten as: 

5<p" + -6<p' = 2X 1 m\5ip, (16) 



8x" + -5x' = (2ml-X 2 m 2 2 )S X 



(17) 
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FIG. 3: The metric function A in the wormhole model. 



FIG. 4: The metric function B in the wormhole model. 



with the exponentially fast damping solutions: 



Sip w C ¥ 



exp ^— -y/2Aimf 
r 

exp ^— yj{2m\ — \2fn\) r 



(18) 



(19) 



where C v , C x are integration constants. Thus the asymptotic solutions go to vacuum ones with the zero energy density 
(Fig. ©)• 



IV. BRANE-LIKE SOLUTION 

In this section we consider a brane-like solution in 4D. Let us chose the metric in the form: 

ds 2 = a 2 (x)(dt 2 - dy 2 - dz 2 ) - dx 2 , 



(20) 



where the metric function a(x) depends only on the coordinate x. This metric describes (2+l)-dimensional spacetime 
embedded in a (3+l)-dimensional spacetime. Using ©, Q and one can obtain the {%) and (*) components of 
Einstein equations (|3|): 



r. 



a \ a / 2 



and the scalar field equations: 



^' + 3-0/ = ^[2 X 2 + A 1 (^-m?)] 7 
a 

x" + 3-x' = x^ + A^-m 2 .)], 
a 

where a prime denotes differentiation with respect to x and the arbitrary constant Vq is choosing as follows 



(21) 
(22) 

(23) 
(24) 

(25) 
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for the purpose of zeroing of a' at x = (see Eqs. (|2T|) and ((26)) below). (Here and further we use the same rescaling 
for all variables and parameters as in the previous section.) 

We will solve the system of equations (|2"T j) - (f24| with the following boundary conditions at x = 0: 

ip(0) - y/3, ip'(0) = 0, 

X (0) = VU, X '(0) = 0, (26) 
o(0) = 1.0, a'(0) = 0. 

The procedure of finding of solutions is the same as in the previous section. The obtained solutions with masses 
77ii ~ 2.59755,TO2 « 3.729 and Ai = 0.1, = 1 for the scalar fields are presented in Fig. (jSJ, for the metric function 
a(x) in Fig. © and for the energy density Tq in Fig. Q. 

We can easily estimate an asymptotic behavior of the solutions. One can see from Eq. 
right hand side tends to zero and the solution of this equation is: 



a = aoe 



that asymptotically the 
(27) 



where ao and a are integration constants. This solution corresponds to the de Sitter-like solution for the space variable 
x. Then, using (|27p and seeking for asymptotic solutions of the equations l|23 p -(|24 p in the form: 



<p = mi- Sip, x = $X, 

where Sip, Sx <C 1 as x — > ±oo, we will have the following equations for Sip and Sx from (|23|) - (f24 

Sip" + iaSip' — 2\xm\8ip 1 

Sx" + 3aS X ' = (2m\-\ 2 ml)5x 



with the damping solutions: 



dip 



C v exp 
C x exp 



-||3a- 



x 



9a 2 + 8Aira 2 



|(3a+ V 9 " 2 + 4 ( 2m i - A 2 TO I) 



(28) 



(29) 
(30) 



(31) 
(32) 



where C V ,C X are integration constants. So we have the solutions that tend asymptotically to the local minimum of 
the potential (|T) at ip = mi and x = 0. 
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FIG. 5: The distribution of the scalar fields <j>, \ near x = 0. 



FIG. 6: The metric function a(x). 
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FIG. 7: The localized energy density Tq(x) near x = with asymptotic anti-de Sitter behavior. 



SPHERICALLY SYMMETRIC SOLUTION 



For consideration of the spherically symmetric problem we take the metric in Schwarzschild coordinates: 



ds 2 



B(r)dt 2 - A(r)dr 2 - r 2 (d6 2 + sin 2 9d(/) 2 ). 



Then the (*) and (£) components of the Einstein equations ([JJ will be: 



-— + — I 

r A 2 r 2 



ir_ 

~B 



1 (B[_ 
2\B 



1 B' 
r~AB 

1 (A_ 
2~A~B 7 



1 

" A 
1 

" A 

&_ 
~B 



= -^ I ^ /2 +x ,2 )-y(^x), 



= 2A 



1 

2A 



and the scalar field equations ([5]) are: 



2B~ 2A 



<p' = A<p[2 X 2 + Xi(<p 2 -mi)] , 



(33) 

(34) 
(35) 
(36) 

(37) 
(38) 



where a prime denotes differentiation with respect to r and in potential (JTJ) Vb = 4 ^2 f° r the purpose of zeroing of 



the energy density as r — > 00. Equation (|36p is consequence of the preceding equations (|34[) - (|35[) . 
Choosing the boundary conditions at r — in the form: 

tp(o) = Vs, p'(o) = 0, 

X (0) = Vo76, X '(0) = 0, 
A(0) = 1.0, B(0) = 1.0 



(39) 



and following the above procedure of obtaining of solutions we find the masses mi « 2.329305, vn,2 ~ 3.0758999 at 
Ai = 0.1, A2 = 1. The results of numerical calculations for the scalar fields are presented in Fig. ([5]), for the metric 
functions A(r),B(r) in Fig. ([9]) and for the energy density Tq in Fig. (fTOjl . 

One can see from (|34 |) -(|35 |) that asymptotic behavior of the metric functions A(r) and B(r) is: 



.4 



1 + a' 



r 



(40) 
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where tq and are constants. Practically ro defines a total mass and B^ the run of time at the infinity. Redefining 
the time variable t we can chose B = 1 as r — > oo. I.e. we have asymptotically flat Minkowski spacetime. The 
corresponding asymptotic scalar field equations (f37|) -(f38 |) with account of 



<p = mi-5<p, x = &X 



(41) 



will be: 



Sip" + -Stp' 
r 

S X " + -Sx' 
r 



with the exponentially fast damping solutions: 



5(p 



2\imf5ip, 

(2m\ ~ ^2ml)5x 



cx:p ( — ^/2Aimf r^j 



cxp (^ — ^(2m\ - A2TO2) r ) 



(42) 
(43) 

(44) 
(45) 



where C v , C x are integration constants. Thus the asymptotic solutions go to vacuum ones with the zero energy density 

(Fig. OH). 

Finally, let us show evolution of the effective equation of state w(r) = p(r)/e(r) where e(r) and p are the effective 
energy density and pressure of the scalar fields. In the case under consideration we have from ([3]) and (|33[) : 



rriO 
J 



e{r) = 
-p(r) 



2A 



(^ 2 + x' 2 )+v^,x) 



2A 



y 2 + x' 2 ) + v(^ x ) 



Then the corresponding effective equation of state will be: 



w(r) — p(r)/e(r) 



^&' 2 + x' 2 ) + v(v,x) ' 



(46) 
(47) 

(48) 



Using the numerical solution obtained earlier, we have the following graph for the equation of state (see Fig. ). 
As one can see from this figure, there is some point r = r* in which the denominator in ()48j) tends to zero (in the 
case under consideration r* ~ 0.956). In the range < r < r* we have w < —1 (w = —1 at r = 0), and in the 
range r* < r < 10 we have w > — 1. At r — > r* from the left w — * —00 and at r — * r* from the right w — > +00. 
Asymptotically w(r) tends to zero, i.e. we have the dust-like equation of state. 

It follows from this that in the entire range < r < r* the WEC is violated, i.e. efficiently the scalar fields describe 
phantom matter. On the other hand in the range r > the WEC is not violated and the scalar fields are non- 
phantom. The scalar fields with the equation of state varying in time are used in the theories describing the present 
acceleration of the Universe (see, e.g., [Hj]). In such theories there is some point of time in which the violation of 
the WEC and transition of usual scalar fields into phantom ones take place. In our case we deal with inhomogeneous 
distribution of the scalar fields with the varying in space equation of state and transition of usual scalar fields into 
phantom ones in some point r = r* . 



VI. CONCLUSION 



We have considered three static solutions for the model of two interacting phantom and ghost scalar fields: the 
model of a traversable wormhole, the brane-like model and the spherically symmetric problem. The self-consistent 
problems with account of the back reaction of the scalar fields on gravitation were investigated. The choice of a 
potential in the form {T]) ensures existence of two local minima that allow to find regular solutions which start and 
finish in one of these minima. In this case the nonlinear problems on evaluation of eigenvalues of the parameters 
mi, m,2 ensuring the mentioned regular solutions were solved. Note that existence of such solutions depends on values 
of the self-coupling constants Ai, A2 and the boundary conditions. There is some range of these parameters appropriate 
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FIG. 8: The scalar fields <f>, \ for the spherically symmetric 



FIG. 9: The metric function A(r),B(r) for the spherically sym- 
metric case. 




w(r) 



FIG. 10: The energy density T$(r) for the spherically symmet- 
ric case. 



FIG. 11: The equation of state w(r) for the spherically sym- 
metric case. 



for existence of regular solutions. Particularly, this range is defining by conditions on existence of the local and global 
minima of the potential (p}: Ai > 0, m\ > \2 r m? i /2\ A2 > 0, m| > \im\/2\ \2tn\ > Ximf. 

For the wormhole model, it was shown that there exist regular solutions in the entire range —00 < r < +00 
with asymptotically flat Minkowski spacetime (see Eqs. I|14p. (fl"5)) ). I.e. the obtained wormhole solution connects 
two flat spacetimes. The radius of the wormhole throat is defined by a minimal value of the function A(r), i.e. 
Rq = min {yl(r)} at r — 0. As one can see in Fig. ([2]), for a remote observer the phantom energy looks like a 

r£( — 00,+Qo) 

compact object localized near the throat of the wormhole with some negative energy density. 

For the brane-like case, the obtained solutions describe (2+l)-dimcnsional spacetime embedded in a (3+1)- 
dimensional spacetime. The asymptotic value of the potential |T]), V(x — ±00) < plays the role of a negative 
cosmological constant, and so the metric (|2"0)) is asymptotically anti-de Sitter with corresponding anti-de Sitter hori- 
zon. In this sense such a solution differs from domain wall solution which has flat asymptotics. 

For the spherically symmetric case we have found particle-like solutions with asymptotically flat Minkowski space- 
time. It was shown that the effective equation of state w(r) for the scalar fields changes essentially along the radius 
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r (see Fig. (JTTJ) ) . There exists some point r — r* dividing the whole space into two regions: in the entire range 
< r < r* , w < — 1 that corresponds to phantom-like behavior of the scalar fields. On the other hand, in the range 
r* < r < oo the equation of state w > — 1 that corresponds to non-phantom (ghost) scalar fields. Similar behavior of 
the equation of state could be obtained both for the wormhole and brane-like models as well. 

Let us note that for all three models, the obtained solutions correspond to soliton-like solutions starting and finishing 
in the same minimum (in the cases under consideration, in the local minimum (f> — mi, x — of the potential l|T|)). 
In models with one scalar field soliton-like solutions exist only in a case of presence of two or more minima and they 
start from one minimum of a potential and tend asymptotically to another one. By the terminology of 17], such 
solutions refer to topological solutions, and our solutions - to non-topological ones. 

Note here that the known constraint on possibility of existence of regular static solutions in spaces with dimen- 
sionality D > 3 does not work in our case. As it was shown in Ref. [3!| (see also Ref. [13)) the mentioned solutions 
for usual (nonphantom) scalar fields do not exist if the potential V > in the whole space. On the other hand, 
for phantom fields the condition V < should be satisfied always (see in this connection Ref. po|). In our case the 
potential |T]) changes its sign that allow to avoid the constraint of the theorem from Ref. [3§| and obtain static regular 
solutions. 

Our attempts to find regular solutions with e = +1 were not succeeded. However, this does not exclude a possibility 
of their existence at some parameters Ai, A2 in the potential ([1]) and some special boundary conditions. 
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